Probing microcavity polariton superfluidity through resonant Rayleigh scattering 
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We investigate the two-dimensional motion of polaritons injected into a planar microcavity by a 
continuous wave optical pump in presence of a static perturbation, e.g. a point defect. By finding the 
stationary solutions of the nonlinear mean- field equations (away from any parametric instability), 
we show how the spectrum of the polariton Bogoliubov-like excitations reflects onto the shape and 
intensity of the resonant Rayleigh scattering emission pattern in both momentum and real space. 
We find a superfluid regime in the sense of the Landau criterion, in which the Rayleigh scattering 
ring in momentum space collapses as well as its normalized intensity. More generally, we show how 
collective excitation spectra having no analog in equilibrium systems can be observed by tuning the 
excitation angle and frequency. Predictions with realistic semiconductor microcavity parameters are 
given. 
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The concept of a quantum fluid has played a cen- 
tral role in many fields of condensed matter and atomic 
physics, ranging from superconductors to Helium flu- 
ids and, more recently, atomic Bose-Einstein conden- 
sates 0. One of the most exciting manifestations of 
quantum behavior is superfluidity, i.e. the possibility of 
the fluid to flow without friction around an impurity 

In this Letter, we investigate the superfluid proper- 
ties of a two-dimensional gas of polaritons in a semicon- 
ductor microcavity in the strong light-matter coupling 
regime 4J. In this system, the normal modes are su- 
perpositions of a cavity photon and a quantum well ex- 
citon. Thanks to their photonic component, polaritons 
can be coherently excited by an incident laser field and 
detected through angularly or spatially resolved optical 
spectroscopy. Thanks to their excitonic component, po- 
laritons have strong binary interactions, which have been 
demonstrated to produce spectacular polariton amplifi- 
cation effects 0, |(j through matter- wave stimulated col- 
lisions Ljjas well as spontaneous parametric instabili- 
ties HUM 

Here, we study the propagation of a polariton fluid in 
presence of static imperfections, which are known to pro- 
duce resonant Rayleigh scattering (RRS) of the exciting 
laser field [HI EE EE E3- We show that superfluidity 
of the polariton fluid manifests itself as a quenching of 
the RRS intensity when the flow velocity imprinted by 
the exciting laser is slower than the sound velocity in 
the polariton fluid. Furthermore, a dramatic reshaping 
of the RRS pattern due to polariton-polariton interac- 
tions can be observed in both momentum and real space 
even at higher flow velocities. Interestingly, the polariton 
field oscillation frequency is not fixed by an equation of 
state relating the chemical potential to the particle den- 
sity, but it can be tuned by the frequency of the exciting 



laser. This opens the possibility of having a spectrum 
of collective excitation around the stationary state which 
has no analog in usual systems close to thermal equilib- 
rium. We show in detail how these peculiar excitation 
spectra can be probed by resonant Rayleigh scattering. 

A commonly used model for describing a planar mi- 
crocavity containing a quantum well with an excitonic 
resonance strongly coupled to a cavity mode is based on 
the Hamiltonian 
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where x is the in-plane spatial position and the field op- 
erators ^x,c( x ) respectively describe excitons (X) and 
cavity photons (C). They satisfy Bose commutation 



rules, [*i(x),*t( x ')] =<S 2 (x- 
tonian h° is: 
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where wcr(k) = u>q \J\ + k 2 /fc 2 is the cavity dispersion 
as a function of the in-plane wavevector k and k z is the 
quantized photon wavevector in the growth direction. 
is the Rabi frequency of the exciton-cavity photon cou- 
pling. A flat exciton dispersion wx(k) = ujx will be 
assumed in the following. In this framework, polaritons 
simply arise as the eigenmodes of the linear Hamiltonian 
Pj l; cj£p(;yp)(k) denotes the dispersion of the lower (up- 
per) polariton branch [FigQJa)]. 
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The external force term proportional to F p describes 
a coherent and monochromatic laser field of frequency 
Lu p (called the pump), which drives the cavity and in- 
jects polaritons. Spatially, it is assumed to have a plane- 
wave profile of wavevector k p — sm9 p u> p /c, 9 P being the 
pump incidence angle, so to generate a polariton fluid 
with a non-zero flow velocity along the cavity plane. The 
nonlinear interaction term is due exciton-exciton colli- 
sional interactions and, as usual, is modelled by a repul- 
sive (g > 0) contact potential. The anharmonic exciton- 
photon coupling has a negligible effect in the regime con- 
sidered in the present study 0. Vx,c( x ) are exter- 
nal potential terms acting on the excitonic and photonic 
fields which can model the presence of disorder. Here, 
results for the specific case of a point defect will be pre- 
sented. Note that point defects can be naturally present 
in state-of-the-art samples ^(J or even be created delib- 
erately by means of lithographic techniques. 



In the quantum fluid language, these are the Gross- 
Pitaevskii equations |2| for our cavity-polariton system. 
For simplicity, an equal rate 7 is assumed for the damp- 
ing of both the excitonic and the photonic fields. In the 
present work, we will be concerned with an excitation 
close to the bottom of the LP dispersion, i.e. the region 
most protected 10] from the exciton reservoir, which may 
be responsible for excitation-induced decoherence [l7j . 

In the homogeneous case {Vx,c = 0), we can look 
for spatially homogeneous stationary states of the sys- 
tem in which the field has the same plane wave structure 
ipx,c( x , t) — exp[i(k p x — u p t)} tpx c as * ne incident laser 
field. The resulting equations 

(wx (kp) -uj p - l - 1+ g 1 rx i 2 ) n + n R r c s = o (4) 

(uc(*p) -Up- ~t) i>c + Mr = -F p , (5) 
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FIG. 1: (a): linear dispersion of the lower (LP) and up- 
per (UP) polariton branches. Typical position of the pump 
wavevector k p (arrow), well in the parabolic region of the LP 
dispersion; the dashed line is the parabolic approximation at 
small k. (b): internal vs. incident intensity curve showing 
bistable behavior. The dashed branch is unstable: the inver- 
sion points A and B are respectively due to a single-mode 
(Kerr) or a multi-mode (parametric) instability. Excitation 
parameters: k v = 0.314/^m -1 , uj p — LULp(k p ) = 0.47meV. 
Cavity parameters are taken from ref . Hal : 7 = 0.1 meV, 
= ua = 1-4 eV, 2 Q R = 5 meV. 



Within the mean-field approximation, the time- 
evolution of the mean fields ipxfii?^ = y&x,c{x)) under 
the Hamiltonian Q is given by: 



are the generalization of the state equation. While the 
oscillation frequency of the condensate wavefunction in 
an isolated gas is equal to the chemical potential /x and 
therefore it is fixed by the equation of state, in the present 
driven-dissipative system it is equal to the frequency to p 
of the driving laser and therefore it is an experimentally 
tunable parameter. As usual, stability of the solutions 
of Eqs. (1 ll'ill has to be checked by linearizing Eq. © 
around the stationary state. For u p > w£,p(k p ), the rela- 
tion between the incident intensity I p cx |-Fj,| 2 and the in- 
ternal one shows the typical S-shape of optical bistability 
(see FigHt and Rfife.jiallillicj). Note that nice hystere- 
sis loops due to polariton bistability have been recently 
experimentally demonstrated m the case k p = 0. In 
the opposite case ui p < o->Lp(k p ) (not shown), the be- 
haviour of the system would instead be the typical one 
of an optical limiter [2J| • 

In the stability region, the response of the system to 
a weak perturbation can be determined by means of the 
linearized theory. In the field of quantum fluids, this ap- 
proach is called Bogoliubov theory By defining the 
slowly varying fields with respect to the pump frequency 
as S(f>i(x-,t) = ^j(x, t) exp(iuipt), the motion equation 
of the four-component displacement vector 5(f>(x.,t) — 
(5<f>x(x-,t), (5(/>c( x , t), 8<j)* x (x, t), 5((>Q (x, t)j reads 
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(3) fd being the source term due to the weak perturbation 
and the spectral operator C being defined as 
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FIG. 2: Left panels: exact Bogoliubov dispersion for the LP 
branches calculated from Eq.Q. Right panels: correspond- 
ing RRS emission pattern in k-space. The intensity has been 
normalized to the transmitted intensity. The k p point (in- 
dicated by a white circle) saturates by far the gray scale. 
Panels (a-f): resonant case A p = respectively in the lin- 
ear regime (a-b), with g\^x\ 2 = 0.2 meV (c-d), lmeV (e-f). 
Panels (g-h): A p > case with g\^x\ 2 = 0.04meV. Panels 
(i-j): A p < case with p|*x| 2 = 0.6 meV. Pump wavevec- 
tor: k p — 0.314 g,m~ l (a-h), 0.408 /an -1 (i-j). Same cavity 
parameters as in FigtTI 

The eigenvalues of the operator C correspond to the 
frequencies of the Bogoliubov eigenmodes of the system. 
For each k, the spectrum is composed of four branches 
u)y P LP (k): for each polariton branch (LP or UP), two 
± branches exist, which are the image of each other un- 
der the simultaneous transformations k — > 2k p — k and 
uj — > 2ojp — uj |2l| . Numerical calculations are shown in 
FigEl For the sake of clarity, only the branches rela- 
tive to the LP have been traced, the ones relative to the 
UP being far away on the scale of the figure. These nu- 
merical results can be understood through the simplified 
analytical approximation that follows. 

Provided the interaction energy g\ipx\ 2 is much 
smaller than the polaritonic splitting uijjp ~^lp, there is 



no significant mixing between the LP and UP branches. 
Since we are interested in nearly-resonant excitation 
close to the bottom of the LP dispersion curve, we can 
describe the system in terms of the LP field i\jlp = 
XlpiPx + C L pipc only, being X LP and C L p the Hop- 
field coefficients quantifying the excitonic and photonic 
components. In the parabolic approximation, cj£p(k) ~ 
ui1 P +hk 2 /2rriLP and the self-coupling constant is g^p = 
g \Xlp\ 4 . The mean-field shift of the polariton mode is 
then Sujmf = ghP IV'l'pI 2 - Under these assumptions, the 
spectrum of the LP Bogoliubov excitations can be ap- 
proximated by the simple expression 



w tP ^ uj p +5k-Vp- — ±^J (2 Sujmf + rjsu - A p )(r)sk - A p ), 

(8) 

where (5k = k — k p , rjsk = fi<5k 2 /2tolp, the flow velocity 
v p = Kkp/rriLp, and the effective pump detuning 

A p = uip - Lu L p(k p ) - Sujmf- (9) 

In the resonant case (A p = 0), the ± branches touch 
at k = kp. The effect of the finite flow velocity v p is to 
tilt the standard Bogoliubov dispersion [2| via the term 
(Jk • Vp. While in the non-interacting case in Fig[^a) 
the dispersion remains parabolic, in the presence of in- 
teractions [Fig'EJc) and (e)] its slope has a discontinuity 
at k = k p : on each side of the corner, the + branch 
starts linearly with group velocities respectively given 
by v l' 1 — c s ± v p , c s being the usual sound velocity of 
the interacting Bose gas c s — ^JHSujmf /ttilp- On the 
hysteresis curve of Fig^b), the condition A p = corre- 
sponds to the inversion point A. If one moves to the right 
of the point A along the upper branch of the hysteresis 
curve, the mean-field shift Sujmf increases and the effec- 
tive pump detuning A p becomes negative. In this case, 
as it is shown in Fig^i), the branches no longer touch 
each other at k p and a full gap between them opens up 
for sufficiently large values of |A p | (not shown). 

On the other hand, the effective pump detuning A p 
is strictly positive on the lower branch of the bistabil- 
ity curve of Fig2Ib). In this case, the argument of the 
square root in JHJ is negative for the wavevectors k such 
that A p > 77,5k > Ap — 2Scomf- In this region, the ± 
branches stick together 15] (i.e. Re[cj + ] = Re[w_]) and 
have an exactly linear dispersion of slope v p (FigEfe). 
The imaginary parts are instead split, with one branch 
being narrowed and the other broadened [pol Elf . For 
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Sujmf > 7/2, that is on the right of point B in Fig^b), 
the multi-mode parametric instability |2l| sets in. In the 
field of quantum fluids, this kind of dynamical instabili- 
ties are generally known as modulational instabilities |'22| . 

The dispersion of the elementary excitations of the sys- 
tem is the starting point for a study of its response to 
an external perturbation. In particular, we shall con- 
sider here a weak and static disorder as described by the 
potential Vc,x( x )- In this case the perturbation source 
term fd is time-independent 

fd = (v x 4% > Vc <t>c , -v x 4% * , -v c 4>c*) T (10) 

as well as the induced perturbation d(j)d = —£~ x ■ fd- 
The static disorder resonantly excites those Bogoliubov 
modes whose frequency is equal to oj p . In the left panels 
of FigEl the excited modes are given by the intersections 
of the mode dispersion with the horizional dotted lines. 
For the specific example of a spatially localized defect 
acting on the photonic component, we have plotted in 
the right panels of FigEl an d in FigOU the photonic in- 
tensity \ipc\ 2 m respectively the momentum and the real 
space for different parameter regimes. These quantities 
correspond to the experimentally accessible far- and near- 
field intensity profiles 0,0,0] of the resonant Rayleigh 
scattering of the pump (i.e. the coherently scattered light 
at the pump frequency lu — oj p ). A similar k-space pat- 
tern is obtained in the presence of a disordered ensemble 
of defects. 
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FIG. 3: Real space resonant Rayleigh scattering emission 
pattern for a localized defect at x — y = acting on the cavity 
photon. The defect potential has a lateral size of 0.8 fim and 
a depth of lmeV. In panel (a), linear regime as in Fig[3Ja) 
and (b). In panel (b), superfluid regime as in Fig|2Ie) and 
(f). In each panel, the emitted intensity normalized to the 
transmitted intensity. 

In the linear regime, the k-space emission pattern 
shown in FigJ^fb) contains a peak at the incident 
wavevector k p , plus the RRS ring [12I fl6| . In the real 



space pattern shown in Fig[3ta), as the polariton fluid 
is moving to the right, the defect induces a propagating 
perturbation with parabolic wavefronts oriented in the 
left direction. 

In the presence of interactions, the RRS circle is trans- 
formed into a oo-like shape with the low-k lobe more 
intense than the high-k one. If A p < 0, the two lobes 
are separated by a gap (Fig|2)), while they touch at k p 
if A p = (Fig|2Jl). In this resonant case, when Suimf is 
large enough for the sound velocity c s in the polaritonic 
fluid to be larger than the flow velocity v p , the slope of the 
+ branch on the low-k side of the corner (Fig|2t) becomes 
positive and there is no intersection with the horizontal 
dotted line any longer. In this regime, RRS is no longer 
possible, and the polaritonic fluid behaves as a superfluid 
in the sense of the Landau criterion |3(. Once normal- 
ized to the incident one, the RRS intensity is strongly 
quenched with respect to the previous cases and no RRS 
ring is any longer present. The weak emission still visible 
in Figl^f) is due to non-resonant processes, which are al- 
lowed by the finite broadening of the polariton modes. As 
no propagating mode is resonantly excited, the perturba- 
tion in real space remains localized around the defect, as 
shown in Fig|3{b). On the other hand, on the bottom of 
the bistability curve (where A p > 0), the polariton gas 
is not superfluid. The RRS intensity is even enhanced 
with respect to the linear regime because of the reduced 
linewidth of the Bogoliubov modes in the regions where 
the ± branches stick together, as shown in Fig[2Kg-h). 

To summarize, the polariton fluid has a superfluid be- 
haviour with respect to static impurities if the equation 
w^p(k) = uj p has no solutions with k ^ k p . If the corre- 
sponding linear regime equation a->Lp(k) = uj p has a set 
of solutions corresponding to the elastic RRS ring, the 
effect of superfluidity is dramatic as the RRS ring is sup- 
pressed. Within the parabolic approximation in Eq. JHJ), 
a simple sufficient condition for superfluidity is found: 
A p < and (rriLpv p — Musmf) < fi\A p \. 

In conclusion, we have shown the strict connection be- 
tween the dispersion of the elementary excitations in a 
quantum fluid of microcavity polaritons and the intensity 
and shape of the resonant Rayleigh scattering on defects. 
In particular, we have pointed out some experimentally 
accessible consequences of polaritonic superfluidity for re- 
alistic microcavity parameters. More in general, thanks 
to the coupling to externally propagating light, micro- 
cavity polaritonic systems appear to be promising can- 
didates for the study of novel effects in low-dimensional 
quantum fluids. 
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